Determination of the dynamical structure of galaxies using optical spectra

Galaxy spectra are a rich source of kinematical information since the shapes of the absorption lines reflect the movement of stars along the line-of-sight. We present a technique to directly build a dynamical model for a galaxy by fitting model spectra, calculated from a dynamical model, to the observed galaxy spectra. Using synthetic spectra from a known galaxy model we demonstrate that this technique indeed recovers the essential dynamical characteristics of the galaxy model. Moreover, the method allows a statistically meaningful error analysis on the resulting dynamical quantities.Comment: 14 pages, 14 figures, Latexfile, MNRAS, in pres

Simulations of the formation and evolution of isolated dwarf galaxies

We present new fully self-consistent models of the formation and evolution of isolated dwarf galaxies. We have used the publicly available N-body/SPH code HYDRA, to which we have added a set of star formation criteria, and prescriptions for chemical enrichment (taking into account contributions from both SNIa and SNII), supernova feedback, and gas cooling. The models follow the evolution of an initially homogeneous gas cloud collapsing in a pre-existing dark-matter halo. These simplified initial conditions are supported by the merger trees of isolated dwarf galaxies extracted from the milli-Millennium Simulation. The star-formation histories of the model galaxies exhibit burst-like behaviour. These bursts are a consequence of the blow-out and subsequent in-fall of gas. The amount of gas that leaves the galaxy for good is found to be small, in absolute numbers, ranging between 3x10^7 Msol and 6x10^7 Msol . For the least massive models, however, this is over 80 per cent of their initial gas mass. The local fluctuations in gas density are strong enough to trigger star-bursts in the massive models, or to inhibit anything more than small residual star formation for the less massive models. Between these star-bursts there can be time intervals of several Gyrs. We have compared model predictions with available data for the relations between luminosity and surface brightness profile, half-light radius, central velocity dispersion, broad band colour (B-V) and metallicity, as well as the location relative to the fundamental plane. The properties of the model dwarf galaxies agree quite well with those of observed dwarf galaxies.Comment: 16 pages, 20 figures, accepted for publication in MNRA

A genetic algorithm for the non-parametric inversion of strong lensing systems

We present a non-parametric technique to infer the projected-mass distribution of a gravitational lens system with multiple strong-lensed images. The technique involves a dynamic grid in the lens plane on which the mass distribution of the lens is approximated by a sum of basis functions, one per grid cell. We used the projected mass densities of Plummer spheres as basis functions. A genetic algorithm then determines the mass distribution of the lens by forcing images of a single source, projected back onto the source plane, to coincide as well as possible. Averaging several tens of solutions removes the random fluctuations that are introduced by the reproduction process of genomes in the genetic algorithm and highlights those features common to all solutions. Given the positions of the images and the redshifts of the sources and the lens, we show that the mass of a gravitational lens can be retrieved with an accuracy of a few percent and that, if the sources sufficiently cover the caustics, the mass distribution of the gravitational lens can also be reliably retrieved. A major advantage of the algorithm is that it makes full use of the information contained in the radial images, unlike methods that minimise the residuals of the lens equation, and is thus able to accurately reconstruct also the inner parts of the lens.Comment: 11 pages, accepted for publication by MNRA

Realistic error estimates on kinematic parameters

Current error estimates on kinematic parameters are based on the assumption that the data points in the spectra follow a Poisson distribution. For realistic data that have undergone several steps in a reduction process, this is generally not the case. Neither is the noise distribution independent in adjacent pixels. Hence, the error estimates on the derived kinematic parameters will (in most cases) be smaller than the real errors. In this paper we propose a method that makes a diagnosis of the characteristics of the observed noise The method also offers the possibility to calculate more realistic error estimates on kinematic parameters. The method was tested on spectroscopic observations of NGC3258. In this particular case, the realistic errors are almost a factor of 2 larger than the errors based on least squares statistics.Comment: 11 pages, 11 figures, accepted for publication by MNRA

Modelling galactic spectra: I - A dynamical model for NGC3258

In this paper we present a method to analyse absorption line spectra of a galaxy designed to determine the stellar dynamics and the stellar populations by a direct fit to the spectra. This paper is the first one to report on the application of the method to data. The modelling results in the knowledge of distribution functions that are sums of basis functions. The practical implementation of the method is discussed and a new type of basis functions is introduced. With this method, a dynamical model for NGC 3258 is constructed. This galaxy can be successfully modelled with a potential containing 30% dark matter within 1r_e with a mass of 1.6x10^11 M_o. The total mass within 2r_e is estimated as 5x10^11 M_o, containing 63% dark matter. The model is isotropic in the centre, is radially anisotropic between 0.2 and 2 kpc (0.88 r_e) and becomes tangentially anisotropic further on. The photometry reveals the presence of a dust disk near the centre

The degeneracy between star-formation parameters in dwarf galaxy simulations and the Mstar-Mhalo relation

We present results based on a set of N-Body/SPH simulations of isolated dwarf galaxies. The simulations take into account star formation, stellar feedback, radiative cooling and metal enrichment. The dark matter halo initially has a cusped profile, but, at least in these simulations, starting from idealised, spherically symmetric initial conditions, a natural conversion to a core is observed due to gas dynamics and stellar feedback. A degeneracy between the efficiency with which the interstellar medium absorbs energy feedback from supernovae and stellar winds on the one hand, and the density threshold for star formation on the other, is found. We performed a parameter survey to determine, with the aid of the observed kinematic and photometric scaling relations, which combinations of these two parameters produce simulated galaxies that are in agreement with the observations. With the implemented physics we are unable to reproduce the relation between the stellar mass and the halo mass as determined by Guo et al. (2010), however we do reproduce the slope of this relation.Comment: Accepted for publication in MNRAS | 12 pages, 8 figure

Wetenschap, op waarde geschat

Oratie uitgesproken door Prof.dr. Sarah de Rijcke bij de aanvaarding van het ambt van hoogleraar in de Wetenschap, Technologie, en Innovatie Studies aan de Universiteit Leiden op vrijdag 17 mei 2019Merit, Expertise and Measuremen